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Abstract 

^v^j ■ In this paper, the densities of the doubly singular beta type I and II distributions 

• ' are found, and the joint densities of their corresponding nonzero eigenvalues are 

. provided. As a consequence, the density function of a singular inverted matricvariate 

' t distribution is obtained. 

o : 

; 1 Introduction 

Let A and B be independent Wishart matrices or, alternatively, let B = YY' where Y has 
a matrix variate normal distribution. Then the matrix variate beta type I distribution is 
defined in the literature as 

f (A + B)- 1 /2 B ( A + B )- 1 /2 j 

U= I B 1 /2(A + B)- 1 B 1 / 2 , (1) 
{ Y'tA + B^Y. 

Analogously, the matrix variate beta type II distribution is defined as 

A _1 / 2 BA~ 1/ ' 2 , 
F= { B 1 / 2 A~ 1 B 1 / 2 ,' (2) 
Y'A _1 Y. 

These can be classified as central, noncentral or doubly noncentral, dependin g on whether A 

and B are central, B is noncentral or A and B are noncentral, respectively; see Diaz-Garcia and Gutierrez- Jaimezl 
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(I2007ll2008bh . addition, such a distribution can be classified as nonsingular, singular or 
doubly singular, when A and B are nonsingular, B is sin gular or A and B are singular, 
respectively, see Dfaz-Garcfa and Gutierrez- Jaimezl ( 2008a ). Under definitions |T]) and © 
and their classifications, mentioned above, the m atrix variate beta type I and II d i stribu - 
tions have b e en stud i ed by different authors; see lOlkin and Rubinl (Il964h . iDickevI (|l967[) . 
Mitral (ll970l) . lKhatril dl970USrivastava and Khatril dl979i).lMuirheadl (11983) . lUhligj (|l994l i. 



Caded (|l996h . iDiaz-Garci'a and Gutierrez- Jaimezl (|l997l . 120061. l2007i l2008allbl ). 



among many 

others. However, the density functions of doubly singular beta type I and II distributions 
have not been studied. In such cases the inverses that appear in definitions ^ and must 
be replaced by the M oore-Penrose inverse. The doubly singular beta type I distribution was 
briefly considered bv lMitra ( 1970h . but this distribution has been the object of considerable 
recent interest. Observe that this distribution appears in a natural way when the number of 
observ ations is smaller than the dimension in multivariate analysis of variance; see lSrivastava 
(|2007h . 

Furthermore, observe that if Z is any random matrix to x r, it is generically called as 
matrix variate if the kernel of its density function is in terms of the trace operator. But 
if the kernel of the density function is in terms of the determinant alone or in terms of 
bot h, then the de term inant and the trace operator, Z, ar e generi cally called matricvariate, 
see IDickevI (|l967l ) and [5 faz-Garcfa and Gutierrez- Jaimezl (|2008ct ). Alternatively, when the 



random matrix A : to x to is symmetric by definition (e.g. Wishart and beta matrices, etc.), 
they are generically referred to as matrix variate. 

In Section [2] of the present study, we review singular matricvariate t distributions and 
determine the distribution of a linear transformation and the joint density function of its 
nonzero singular values. In Section [3l an expression is provided for the density function of 
the doubly singular matrix variate beta type II distribution, and the joint density function 
of its nonzero eigenvalues is determined. 

Similar results are provided for the doubly singular matrix variate beta type I distribu- 
tion; see Section |4j Finally, in Section [5l we study the inverted matricvariate t distribution, 
also termed the matricvariate Pearson type II distribution. 



2 Preliminary results 

Let C r ^ m {q) be the linear space of all to x r real matrices of rank q < min(TO,r) and let 
C^ m (q) be the linear space of all to x r real matrices of rank q < min(TO, r), with q distinct 
singular values. The set of matrices Hi E C r ,m such that H^Hi = is a manifold denoted 
as V rym , termed the Sticfcl manifold. In particular, V m . m is the group of orthogonal matrices 
0(m). The invariant measure on a Stiefel manifold is given by the differential form 

r m 

(HidH!) = /\ f\ h'jdht 

2—1 j — 

written in terms of the exterior product (A), where we choose an m x (to — r) matrix H2 

such that H is an m x m orthogona l matr ix, with H = (Hi:H2) and where dh is an m x 1 
vector of differentials; see iMuirheadl (| 19821 . Section 2.1.4). Moreover 

r or —mr/2 r <ym m 2 /2 

/ (H[dH 1 ) = - 1 — and / (H'dH) = -— — - (3) 

iH lS v,, m r r [TO/2J i He o(m) r m [m/2J 

Denote by <S m , the homogeneous space of to x to positive definite symmetric matrices; 
and by S m (q), the (mq — q(q — l)/2)-dimensional manifold of rank q positive scmidefinitc 
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to x to symmetric matrices; and by <S+(g), the (mq — q(q — l)/2)-dimensional manifold of 
rank q positive semidefinite to x m symmetric matrices with q distinct positive eigenvalues. 
Assume A G S m ; then chj(A) denotes the i-th eigenvalue of the matrix A. Moreover, let 
A + and A - be the Moore -Penrose inverse and any symmetric generalized inverse of A, 
respectively; see [Rao (1973). Finally, A 1 / 2 is termed a non-negative definite square root of 
A if A 1 / 2 is such that A = A 1 / 2 A 1 / 2 . 

Using the notation of iDiaz-Garcia et al. I (| 19971 ) fo r the singular matrix variate nor- 
mal an d Wishart and pseudo-Wishart distributions, from lDi'az-Garcfa and Gutierrez- Jaimez 
(|2008d) we have: 

Lemma 2.1. Let T G £+ m (r s ) be the random matrix 



A 1 / 2 ) Y + v 



where A 1 / 2 A 1 / 2 = A G 5+(q) ~ VU* (n, 0), © G S m (r @ ) with r & < to and q = min(TO, n), 
independent of Y G £tm{ r s) ~ ^mxr (fli ^™ ® " e ^v(V s ) wii/i r s < r and m > q > 
r H > 0. TTien </ie singular matricvariate T /ias i/ie density 



i=i 
with 



c(m,n,q, gi , r e , r H , r J 
ch^S)" 1 ^ Jlch^©)"/ 2 



fjchz [0- + (T-m)S-(T-/x)'] 



,-,-(n+r s )/2 



(dT), 



(4) 



2=1 



c(m,n,q,q 1 ,r & ,r 3 ,r a ) 



7r n(9-f e )/2-(n-l-r a )(«i-r c ,)/2-mr s /2p^ J( n + r a )/2] 
2(mr H +nr e )/2-(n+r s )r Q /2p^ 



(5) 



where q\ = min(TO,n + r B ), r a = rank[& + (T — fi)S (T — fi)'] < to, and (dT!) denotes 
the Hausdorff measure. 

In particular observe that if S = I r , i.e. r H = r, r@ = r a = m, & = I m , q = n and 
<7i = n + r, we have 

dF T (T) = ^-K^ + 2n)/2 rn + .[(n + r)/2] |im + (T _ ^ )(T _ ^y,- (»-Hr)/ 2(dT) (6) 
where T G £+ m (r) and now (dT) denotes the Lebesgue measure. 

Theorem 2.1. Under the condition of Lemma \2.1[ let [i = and X = TC + G £+_ m (r H ), 
where B = CC' swc/i t/ia£ C G £+_ r (r B ). TTien i/ie density function o/X is 



c(m,n,q, q 1 ,r @ ,r s ,rj -q ^ 

nch J (©)«/2 



XX 



A-(n+r s )/2 



(dX), 



(7) 



where q± = min(m, n+r H ), r a = rank(& + XX') < to, and now (dX) denotes the Lebesgue 
measure. 

Proof. Make the change of variables T = XC as described by IDiaz-Garcia ( 2007 ) 

(dT) = JJchi(CC') ro/2 (dX) = nch l (H)™/ 2 (dX). (8) 



i=l 
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Also, observe that 

C+HC+' = C+(CC')C + ' = C+C(C+C)' = C+C = I rs , (9) 

since C and C + C have the same rank r s and C + C is of the order r H x r s . Finally, note 
that by Q 

TS-T' = XC'S-CX' = X(CTEC'-)-X' = X(C+HC+')-X' = XX'. (10) 
By substituting ((8|) and (fTOjl in (J4j) we obtain the desired result. □ 

Now let Ki, i = 1, . . . , r be the nonzero singular value of T. Then taking /i = in ^ we 
have 

Theorem 2.2. The joint density function of K\, . . . , n r , the singular values of T is 



r(2n—m) /2 

2 r r n+r [(n + r)/2] 



r B [»/ 2 ]r r [r/2]r r R^n (1 g ^ «*)[/\ d *j (") 

where n\ > ■ ■ ■ > n r > 0. 

Proof. Let T = QiDtP' be the nonsingular part of the singular value decomposition 
(SV D), where Qi <£ V r .m , P € and D T = diag(«i, . . . , /%■), K\ > ■ ■ ■ > n r > 0. Then, 



from iDiaz-Garcia et al. ( 19971 ) 



(dT) - 2- r |D T | m ^n (k? - (dD T )(QWQ 1 )(P'dP). 
The result is then obtained immediately from ([S]), noting that 

\lrn + QiDtP'PDtQ'!! = |I, + T> 2 T \ 

and integrating over Qi e V r ^ m and P S O(r) using ([3]). □ 



3 Doubly singular beta type II distribution 

Theorem 3.1. Let A 1 / 2 A 1 / 2 = A e S+(<j A ) ~ W™ A (n, 0), £ 5 m (r e ) luii/i r e < m 
and g A = min(m, n), independent of Y e £^ m (7" s ) ~ A/" m x r s (0,I m ® S), S € 5r(f* s ) 
r s < r and m > q > r H > 0. The doubly singular matrix variate beta type II distribution is 
defined as 

F = (A 1 / 2 )+(YH-Y)(A 1 / 2 )+ = (A 1 / 2 )+B(A 1 / 2 )+ 
where B = (YE~Y) 6 Sm( r =) ~ 'PWrn (r B ,lm)- The density function o/F is 

c(m, n, g, ft , r,, r. , r J | ^ | (r . g _ m _ 1)/2 ^ ^ (@ _ + ^-(n+r.J/a ^ (12) 

i=i 

where F = QiDfQi is i/ie nonsingular part of the spectral decomposition, with Qi G V r= ,m 
and Dp = diag((5i, . . . , 8 r _), 8\ > ■■■ > 5 r _ > and where (dF) denotes the Hausdorff 
measure. 
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Proof. By (fit))) , observe that 

F = (A 1 / 2 )+(YH-Y)(A 1 / 2 )+ = TETT' = XX' 

Let X = QiDxP' be the nonsingular part of the singular value decomposition (SVD), where 
Qi € V r m , P £ 0(r) and D x = diag(«i, ...,«,.), «i >■••> «v > 0. Then F = XX' = 
QiD xP'PDxQi = QiDyQi = QiDpQi, with D F = diag(<5i~ . . , S r ) = D 2 -. Then, 



from Dfaz-Garcia et al. 



\uam, 



(dX) = 2- r ^\T> F \^~ m ^/ 2 (dF)(P'dP). (13) 

Substituting ([T3"j) in J7]) and integrating over P £ 0(r) using Q gives the stated density 
function for F. □ 



In particular, if © = I m , q\ = n + r H and q = n, then r a = m and 

r r 'T B /!)lr + r B) /91 ] V F \ ir -- m - 1)/2 \lrn + F r ("+^)/ 2 (dF). (14) 

An alternative way to obtain (| 14[) when S = I r and therefore r s = r is by using the 
third definition of F given in This denotes the matrix variate as F. 

Theorem 3.2. Let A 1 / 2 A 1 / 2 = A £ S+{n) ~ PW^ {n,l m ), independent ofY £ C+ m {r) ~ 
A/" mX r(0,I m ® I r ), and m > n > r > 0. T/ie doubly singular matrix variate beta type II 
distribution is defined as 

F = Y'A+Y, 

which has the density function 

^-r(r+2n- m )/2 r n+r[(n + r)/2] |p| (m-r- 1)/2 |^ + p | - (n+r)/2 f^p\ ^ 

r„[n/2]r r [m/2] 

where (dF) denotes the Lebesgue measure. 
Proof. Observe that 

F = Y'A+Y = T'T, 

where the density function of T is given by ([6]). Now as in the proof of Theorem 12. 2i let 
T = QtDtP'; hence F = T'T = PD T QiQiD T P' = P'D f P', with D f = D 2 ,. Then by 



Diaz-Garcia et al. I (|l997l ) 

(dT) = 2- r \F\ m - r - 1 ^ 2 (dF)(Q' 1 dQ 1 ). (16) 

The proof is completed by substituting (fTB)) in © and integrating over Qi £ V r , m using ©. 
□ 

Note that T heorems 13.11 and 13.21 are generalizations of Theorems 10.4.1 and 10.4.4 in 



Muirheadl (|1982T ) for the central doubly singular case. 

Theorem 3.3. Let 5i,...,S r be the nonzero eigenvalues ofF (with density function \1J$ 
and let r s = r) or equivalent be the nonzero eigenvalues of F (with density function U5\) ). 
Then the joint density function of Si , . . . , S r is 

^ ri^^ n n « - « (A *) <"> 

where 8\ > ■ ■ ■ > S r > 0. 
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Proof. This result can be obtained by any of the following methods: 

i) By making the transformation F = QiDpQi in (| 14[) with r s — r and integrating over 

Qi & V„„ using ([3]). 

ii) Alternatively, by making the transformation F = PDpP' in (TTSJ) and integrating over 

P e O(r) using ([3]), the proof is completed. 

1 /2 

iii) A third proof is derived immediately from Theorem 12.21 observing that Ki = 5 i , i = 

1 , . . . , r and so 

/ r \ r 

-1/2 

\i— 1 / i— 1 \ «— 1 / 

Remark 3.1. Finally, from Mitra ( 197(1 Lemma 3.1) observe that 



f\d Kl )=2- r l[s; 1/2 \f\dsA □ 



a) The density function (| 14[) is invariant under any arbitrary matrix G S m and Y G 

b) Similarly, Theorem 13.21 is invariant if A G *S+(n) ~ WJJ, (n, 0), independent of Y e 

£+mM ~ Knxr(0, © ® Ir), where G S m . 

c) Therefore, Theorem 13.31 too, is invariant under conditions established in items a) and 

b). 

4 Doubly singular beta type I distribution 

Theorem 4.1. Let A 1 / 2 A 1 / 2 = A G S+(n) ~ PW™ (n,I m ) &e independent of Y G 
£+ m (r s ) ~ A/" TO x^(0,I m (8> 3), 2 G S r (r s ) with r s < r and m > q > r s > 0. The 
doubly singular matrix variate beta type I distribution is defined as 

U = ((A + YE-Y) 1 / 2 ^ (YETY) (a + YETY) 1 / 2 

= fA + B) 1 / 2 ) + BfA + B) 1 / 2 ^ + 



where B = (YH Y) G Sm( r s) ~ T^Wm (r H ,I m ). TTie density function o/U is 
r„ +r= [(n + r s )/2] 



r„[n/2]r r= [r s /2] 



| Du |(r s -m-l)/2|j m _ V \(n-m-l)/2^ dX j^ ^gj 



where U = QiDuQi *s i/ie nonsingular part of the spectral decomposition, with Qi G V r _ !m 
and Du = diag(Ai, . . . , A r _ ), 1 > Ai > • • • > A r _ > and where (dXJ) denotes the Hausdorff 
measure. 

Proof. From lPiaz-Garcia and Gutierrez- Jaime d ( 20061 ) it is known that if F G 5+ (r H ) is a 
singular matrix variate beta type II distribution, then the matrix variate U = I m — (I m +F) _1 
has a singular beta type I distribution; moreover 

(dF) = |I rs - Dur(" l+1 -' S )/ 2 |I m - U|-( m+1+r ^/ 2 (dU) 

where U = QiDuQi with Qi G V r=im andDu = diag(Ai, . . . A r= ), 1 > Ai > • • • > A r > 0. 
The proof follows from making the transformation F = (l m — U) -1 — \ m in (fT4|). noting 



that F = (I m — U) 1 U and Dp = — Du) 'Du; see lDiaz-Garcia and Gutierrez- Jaimez 



(|2006h . = □ 



G 



Theorem 4.2. Let A 1 / 2 A 1 / 2 = A 6 S+{n) ~ 7>W™ (ra,I m ), independent o/Y € £+ m (r) ~ 
A/ r mX r(0,I m ® I r ), and m > n > r > 0. The doubly singular matrix variate beta type II 
distribution is defined as 

U = Y'(A + YY') + Y, 

and its density function is given by 

^-r(r+2n-rn)/2 r n+r[{n + r)/2] I I (m-r- 1) /2 jj^ _ U T ( ™ -m ~ 1)/2 (rfU) , (19) 

r„[n/2]r r [m/2] 



where (dU) denotes the Lebesgu 



e measure. 



Proof. The proof is obtained immediately by making the transformation F = (I m — 
U)- 1 !! with (dF) = |I m -U|-( r+1 ) in (TT5|). □ 



Theorem 4.3. Let U € 5m( r ) ^ e a matrix variate with density function H8\) and r s = r 
(or let U £ iS r , matrix variate with density function h!9fl ). Then the joint density function 
of X\, . . . ,X r , the nonzero eigenvalues o/U (or\5), is 

r (2n-m)/2 r» +r [(n + r)/2] 

r„[n/2]r r [r/2]r r [m/2] 

f[ A| m -- 1)/2 (1 - X^-^HiX, - X 3 ) (A dA,) (20) 

2 = 1 i<j / 

where 1 > Ai > • • • > A r > 0. 

Proof. The proof can be obtained in any of the following ways. 

i) By making the transformation U = HiDuH' in (|18[) with r s = r and integrating over 

Hi G V r ,m using (0J). 

ii) Alternatively, by making the transformation U = GDgG' in (| 19|) and integrating over 

GeO(r) using ©. 

iii) A third proof is derived immediately from Theorem l3.3l observing that Ai = (1 — ) 1 (5.^ , 

i = 1, . . . , r, and so 



Finally, observe that conclusions analogous to those set out in note 1 can be obtained 
for the case of the matrix variate beta type I distribution. 



5 Matricvariate inverted t distribution 

Let us now find the matricvariate inverted t distribution, also termed the matricvariate 
Pearson type II distribution. In this case, we determine its density function from the doubly 
singular beta type I distribution (fl8|) . 
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Theorem 5.1. Let R € Cj m {r) be the random matrix, 



R 



(A + YY 



1/2 



w/iere A 1 / 2 A 1 / 2 = Ae5+(n)~ ?W m (n, I m ), independent ofY G £+ m (r) ~ 7V„ IXI .(0. I m ® 
I r ) and < r < n < m. TTien i/ie matricvariate R /ias t/ie density function 

dF R (R) = ^-r(r + 2n)/ 2 rn + ,[(n + r)/2] |lm _ rr/| _ („ +r ) /2(dR)) (21) 
where I m — RR' > and (dR) denotes the Lebesgue measure. 



Proof. Observe that U = RR' = (A + YY') + j] 1//2 YY' (A + YY')^] 1 / 2 . Then 
from IDiaz-Garcia et al. (1997), 



(dR) = 2~ r |Du 



(j--m-l)/2 



(dU)(G'dG), 



where R = HiD R G' and U = RR' = HiDrG'GDrH; = HiD^H^ = HiDuH'j with 
Du = D^, Hi e V r ,m, G G 0(r) and D R = diag(ri, . . . r r ), 1 > t x > ■ ■ ■ > r r > 0. Thus 

(dU) = 2 r |Nr( r - m - 1 )/ 2 (dR)(G'dG)- 1 . 
On substituting in (fT5|) and integrating over G € 0(r) using ©, the proof is complete. □ 



Alternatively, Theorem 15.11 can be proved from the distribution of U in an analogous 



way. 



Conclusions 



IDiaz-Garcia and Gutierrez- Jaimez ( 2008c ) studied the density of the singular matricvari- 
ate t distribution and its application to sensitivity analysis. In this work we study some 
properties of this distribution, namely the distribution of a linear transformation and the 
j oint density function of its singular values. Uhlid (1994), Diaz-Garcia and Gutierrez- Jaimez 
(| 19971 . l2008al ) studied the singular matrix variate beta type I and II distributions in central 
and noncentral cases. Now we obtain the doubly singular matrix variate beta type I and 
II distributions in the central case and find the joint density function of their eigenvalues. 
Finally, we obtain the central density function of the Pearson type II matricvariate. 
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